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ABSTRACT 

Studies of outer regions of spirals disks are fundamental to our understanding of both the process of galaxy 
assembly and the subsequent secular evolution of galaxies. In an earlier series of papers we explored the extent 
and abundance gradient in the outer disk of NGC 300 and found an extended purely exponential disk with a 
metallicity gradient which flattens off in the outermost regions. We now continue the study of outskirts of pure 
disk spirals with another Sculptor Group spiral, NGC 7793. Using Gemini Multi Object Spectrograph camera 
at Gemini South, we trace the disk of NGC 7793 with star counts out to ^ 9 scale lengths, corresponding to 
1 1.5 kpc at our calibrated distance of 3.61 ± 0.53 Mpc. The outer disk of NGC 7793 shows no evidence of a 
break in its light profile down to an effective surface brightness of ^ 30 mag arcsec"^ 3 mag arcsec"^ deeper 
than what has been achieved with surface photometry) and exhibits a non-negative abundance gradient within 
the radial extent of our data. 

Subject headings: galaxies: abundances — galaxies: individual (NGC 7793) — galaxies: stellar content — 
galaxies: structure 



1. INTRODUCTION 

The study of the extreme outer disks of nearby spirals has 
produced some remarkable findings in recent years. These 
advances have been brought together by a confluence of ob- 
servational and theoretical discoveries. On the observational 
side, star counts offer a superior method for tracing faint stel- 
lar populations in the outskirts of galaxies, compared to tra- 
ditional surface photometry ([ Pritchet & van den Bergh ll 19941: 
Ferguson et al. 2002; Bland-Hawthorn et al. 2005; Irwi n et al] 
2005; de Jong et al. 2007; Radburn-Smith et al. 2011). While 
surface photometry is limited by diffuse background sources 
at the level of ~ 27 mag arcsec"^, observations of re- 
solved stellar populations do not suffer from the same lim- 
itations and allow us to reach 3-5 mag arcsec"^ deeper 
in effective surface brightness. On the other hand, the N- 
body/SPH and semi-analytical models of disk formation are 
now able to address the importance of secular effects in the 
evolution of spirals. Sc attering of stars by t ransie nt spi- 
ral arms first studied by Sel lwood & BinnevI (l2002h could 
potentially have a significant effect on evolutio n of spi- 
ral disks (Roskaretal. 2008b a; Schonrich & Binnev 2009; 
fSdnchez-Blazquez et al. 2009; Mar tinez-Serrano et al.. 2009; 
iMinchev & Famaey 2010; Minch ev et all 120 lOl) . The ef- 
fects of the secular evolution of galaxies are expecte d to be 
most prominent in the outermost parts of disks (Rosk aret alj 
l2008a). Additionally, due to their long dynamical timescales, 
the outer regions of spirals have largely retained fossil record 
from the epoch of galaxy assembly in the form of spatial dis- 
tributions, kinematics, ages and metallicities of their stars. 
These factors make outskirts of spirals particularly useful for 
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testing models of galaxy formation and evolution. 

Predictions that models of secular evolution o f spirals make 
for stellar ages and metallicities in disk o utskirts (iRoskar et alJ 
I2008blal : ISdnchez-Blazquez et alJ 120091) are still to be con- 
firmed in a sufficiently large sample of galaxies. Age be- 
havior is very difficult to observe in galaxies too dis tant for 
the fu ll sta r formation history to be modeled. Barker et alj 
(120071) and I Williams et al] (l2009h demonstrate nicely how an 
age gradient can be derived over the whole disk of M33; how- 
ever, this kind of work is only possible for the few closest 
spirals. On the other hand, abundance gradients in galax- 
ies can be studied relatively easily out to a few Mpc using 
broad-band photometry. It has been well-established that in a 
disk which is growing inside-out, metallicity decreases from 
the center of the galaxy o utwai'd (iGoetz & Koeppenlll992t 
iMatteucci & Francoislll989h . However, an increasing num- 
ber of galaxies show signs of an abundance gradient flattening 
in their outermost regions ( Andrievsky et al. 2004; Yon g et al. 
|2006; Carraro et al. 2007; Pedicelli et al. 2009; Bres olin et al. 
|2009b; Vlajic et al. 2009). WTiile models mentioned above 
offer possible explanations for the flattening, its exact cause 
is not yet clear and abundance data on more galaxies of differ- 
ent sizes and outer disk structure will help explain the slope 
of abundance gradient in the outermost regions of spirals. 

As discussed in our earlier papers, the Sculptor Group is an 
ideal laboratory for studies of pure spirals: (i) its high galac- 
tic latitude minimizes foreground extinction and contamina- 
tion from Milky Way stars, (ii) most of its galaxies are iso- 
lated disk systems, and (iii) the distance of the Sculptor Group 
(2-4 Mpc) allows for the resolved stellar population studies 
from the ground. At the far side of the Sculptor Group (^ 3.6 
Mpc), NGC 7793 is a late-type spiral with an absolute mag- 
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Fig . 1 .— DSS wide-field image of NGC 7793 with GMOS fields overlaid. 

nitude M = -18.46, slighdy brighter than NGC 300 (-18.12; 
absolute magnitudes were derived usin g corrected total ap- 
parent magnitude from Carignan' ("1985) and the distance de- 
rived in ^3.1b . Carignan (1985) studied the light profile of 
NGC 7793 out to 6' using surface photometry and derived 
a scale length of 1.10^ ^^4^ ^Pc) in B. A neutral hydrogen 
study of lCarignan & PuchS (119901) reveals an H I disk extend- 
ing out to 11.2' (11.8 kpc), with a surface density at the last 
data point of 0.01 M© pc"^. 

The plan of the paper is as follows. In ^ we describe 
the observations, data reduction, photometry and complete- 
ness analysis. Next, we present results for the distance of 
NGC 779 3 as estimated from the tip of the red gia nt b ranch 
(TRGB, ^TTT i. color-magnitude diagram (CMD, ^12] |, star 
counts profiles ( ^3.3t, sur face brightness profile ( ^3.4t and 
metallicity gradient 03.5b . Discussion and conclusions fol- 
low in H 

2. OBSERVATIONS, DATA REDUCTION AND PHOTOMETRY 
2.1. Observations and Data Reduction 

The data were obtained using the Gemini Multi Object 
Spectrograph (GMOS) on the Gemini South telescope over 
three nights in 2005 August as a part of the program GS- 
2005B-Q-4. Deep g' and /' images of two major axis fields 
on each side of the galaxy were taken; the locations of the 
fields (SE and NW) are shown in Figure [T] GMOS field of 
view is 5'.5 on a side. (At the distance of NGC 7793 (3.61 
Mpc, 93.11 1. 1' corresponds to 1.05 kpc.) The average FWHM 
of the data is 0".6 (SE) and 0".8 (NW) in g', and 0".5 (SE) 
and 0".7 (NW) in/' band. 

To reduce the data we employed the standard IRAF/Gemini 
routines which included (i) bias subtraction and flat fielding 
(gi reduce), (ii) (for /' data only) creating the master fringe 
frame from the individual reduced frames (gi fringe) and 
subtracting it from the individual images (girmf ringe), 
(iii) mosaicking of individual GMOS CCDs into a single ref- 
erence frame (gmosaic), and (iv) combining the dithered 
exposures into a final image (imcoadd). 

We obtained 1 3 x 600 s exposures in g' and 22 x 600 s ex- 
posures in /' band per field, bringing total on-source exposure 
time to 1 1 .7 hours. The data for the SE field were taken during 



the nights of 2005 Aug 09UT (hereafter: first night) and 2005 
Aug lOUT (hereafter: second night). Only 5 (out of 13) g'- 
band and 21 (out of 22) /'-band images were observed on the 
first night and the observing log indicated that a thin cloud 
might have affected /'-band observations. In addition, only 
one standard star field observation was recorded. Remain- 
ing science images of SE field, as well as three standard stars 
fields were observed the following night. The bulk of the data 
of the NW field was coflected the night of 2005 Aug 1 lUT (12 
g'-band and 21 /'-band images; hereafter: third night). How- 
ever, no photometric standard stars were observed that night. 
The remaining science frames (one in each of the bands) were 
observed the previous night under photometric conditions. 

Initial analysis of the photometry revealed a suspicious dis- 
crepancy in /'-band magnitude distribution between the two 
fields. In correspondence with the Gemini staff it was con- 
firmed that this was most likely due to the non-photometric 
conditions on the night of 2005 Aug 09UT, when the majority 
of the science frames of the SE field were taken. Accounting 
for this, and the fact that the photometric standard stars ob- 
servation were only taken on the second night of the run, we 
decided to proceed in the following manner: 

- SE field, /'-band: combine 21 images taken on the first 
night and compare photometry of this deep image with the 
photometry extracted from the single image observed on the 
second night. 

- SE field, g'-band: combine separately 5 images observed 
on the first, and 8 images taken on the second night and com- 
pare the photometry between the two images. 

- NW field, /'-band: compare the photometry extracted 
from the deep image created by combining 21 images taken 
on the third night of the run with the single image taken on 
the second night. 

- NW field, g'-band: similarly to the /'-band case, compare 
a combined deep image created from the 12 exposures taken 
on the third night, with the single image taken on the second 
night of the run. 

In the following sections, we will refer to the (combinations 
of) images taken on the second night of the run as calibration 
images and to the ones observed on the first and third night as 
final images. 

2.2. Photometry and Artificial Stars Tests 

To extract stellar photometry we used the standalone ver- 
sion of DAOPHOT and ALLSTAR packages (Stetson 1987). 
Following the initial runs of FIND and PHOTOMETRY routines, 
which were used to catalog objects in the image and mea- 
sure their aperture photometry, we proceeded to determine 
the point spread function (PSF) for each image. Depending 
on the filter and the field, 80-220 moderately bright isolated 
stars in each field were selected as PSF stars and used to iter- 
atively compute the PSF. The PSF stars were "hand-picked" 
and their radial, contour and mesh profiles were visually ex- 
amined within the IRAF/DAOPHOT package. The calculated 
point spread function was used to subtract the PSF stars from 
the original image; the positions of the subtracted PSF stars 
were inspected again and the PSF stars which did not subtract 
cleanly were excluded from the PSF calculation. In addition, 
stars with subtraction errors which differed more than 3 a from 
the mean value were also excluded. The next iteration of the 
point spread function was calculated using images in which, 
within the fitting radius of each PSF star, all but PSF stars 
have been subtracted. This was followed with yet another vi- 
sual inspection within IRAF/DAOPHOT as described above. 



The Outer Disk of NGC 7793 



3 



The whole procedure was repeated once more to derive the fi- 
nal PSF. Finally, ALLSTAR was used to fit the calculated PSF 
to all stars in the object catalogs and determine their photom- 
etry. 

Comparison of the final and calibration images revealed 
the following. In the /'-band, we found the difference in 
photometry extracted from the calibrated and final images of 
0.396 ±0.012 for the SB field and 0.073 ± 0.006 for the NW 
field. In the g'-band we found, within uncertainties (2a), no 
difference in magnitudes between calibration and final im- 
ages. We therefore decided, in the subsequent analysis, to 
use the master frames created by combining the full sample 
of 13 frames in g', and a reduced sample of 21 images in the 
/'-band. Two /'-band images, one of each field, taken on the 
second night were used to calibrate the photometry of the final 
images which were very likely taken under non-photometric 
conditions. 

We performed photometric calibration using the standard 
stars observations taken during the second night of the observ- 
ing run. The derived zero points were in excellent agreement 
with those stated on the Gemini website (28.31 (our value) 
vs. 28.33 (GMOS value) for g' and 27.92 (our value) vs. 
27.93 (GMOS value) for /' band). In addition, we applied 
additional offset (0.396 for the SE field, 0.073 for the NW 
field) to the photometry derived from our final /' images, as 
explained above. 

The DAOGROW software suite (IStetsonll 19901) was used to 
calculate aperture corrections. We subtracted from each im- 
age all but PSF stars and determined their aperture photome- 
try in a series of apertures. DAOGROW uses the information 
on aperture photometry to derive the "total" magnitude of a 
star The weighted mean difference between the PSF-based 
magnitude (from ALLSTAR) and the "total" magnitude (from 
DAOGROW) of the PSF stars was adopted as 'aperture cor- 
rection' and applied to PSF-based magnitudes of all stars in 
the frame. The uncertainties in the aperture correction were 
in the range 0.006-0.008. 

Completeness, crowding and photometric uncertainties of 
the data were assessed using artificial star tests, described in 
more detail in Vlajic et al. (2009). We added 900 artificial 
stars to each frame in a 30 x 30 grid with the cell size of 
69 pixels (10"); for each combination of field and filter 100 
artificial star test runs were performed with the grid origin 
randomly offset between the runs. Photometry (i.e. magni- 
tude range and color distribution) of artificial stars was de- 
termined by randomly sampling the original color magnitude 
diagrams. Artificial stars were added to the frames using 
the DAOPHOT/ADDSTAR routine and the resulting images 
were analyzed with the identical data reduction pipeline as 
the original frames. We consider an artificial star recovered 
if it is detected in both filters with the difference between in- 
put and recovered magnitude < 0.5 mag. Completeness as 
function of magnitude is calculated as a ratio between the 
number of recovered and input stars in a given magnitude 
bin. Our photometry is 50% complete down to (§',/')=(27.35- 
27.47,26.75-27.10). We also calculate radial completeness as 
a ratio of the number of recovered and input stars in a given 
radial bin; the information is later used to correct radial star 
counts and effective surface brightness profiles for crowding 
in the innermost regions. Furthermore, we used information 
on the CHI and SHARP parameters of recovered artificial star 
tests to discard spurious objects from the ALLSTAR catalog. 
CHI and SHARP are often employed as indicators of the non- 
stellar nature of an object; outHers in the plots of CHI and 



SHARP as a function of magnitude can be relatively safely as- 
sumed to be semi-resolved galaxies, blends, cosmic rays or 
image defects. We removed from further analysis all stars ly- 
ing outside the envelope demarcating 3 standard deviations 
from the artificial stars' mean CHI and SHARP values. These 
cuts, together with a cut on error in color (w 0.15) remove 
^ 50% of objects from the original matched catalogs. 

3. RESULTS 

3.1. Tip of the Red Giant Branch Distance 

Stellar evolution models (e.g. ISalaris & Cassisil 119971: 
Iben & Renzini 1983) show that below a metallicity of 
[Fe/H]« -0.7 the brightness of the tip of the red giant 
branch (TRGB) is expected to be roughly constant, re- 
gardless of stellar age and m etallicit y. Observat i onally , 
iDa Costa & Armandroffl (Il990t) and iBellazzini et al.1 (1200 lb 
confirmed that for a sample of Galactic globular clusters span- 
ning a large range of abundances (-2.1 < [Fe/H] < -0.7) 
and ages (2-15 Gyr), the absolute /-band magnitude of the 
TRGB is stable and insensitive to age and metallicity, which 
enables its application as an extragalactic distance estimator. 
In general, there is been a good agreement between the dis- 
tances obtained using the TRGB and classical methods such 
as Cepheids (Sakai et al. 1996; FerrtireseeJ^al. 2000). 

Following lLee et al.l(il993h and lSakai et'al.l (119961) . we em- 
ploy our deep /'-band photometry of two fields in the outskirts 
of NGC 7793 to determine the galaxy's distance. Figure |2] 
shows the result of this exercise. In order to increase the sam- 
ple of stars used in the detection, we combined photometry of 
the two fields. This results in a value consistent with the one 
we derive from each field separately, while having a narrower 
and better defined peak. 

The top panel of the Figure shows the /' band luminosity 
function, with (solid line) and without (dashed line) the com- 
pleteness correction. In correcting for incompleteness, we 
only consider stars that are above 50% completeness limit in 
both g' and /'. In the middle panel, the smoothed luminosity 
function is shown. In calculating this function, discrete stellar 
magnitudes are replaced by Gaussians with the width equal 
to the star's photometric error, according to the expression 
dS^aTet al.. l996 ): 



1 

$(m) = V exp 



i=l 



(nii — m) 



2a 



A smooth adaptive edge-detection filter of the form: 
E(m) = $(;« + CT„,) - $(m - a„) 

is then applied. Here, ct,„ is the mean photometric error of 
all stars in the magnitude range [m-0.05,m + 0.05], as de- 
termined from artificial star tests. The maximum of the out- 
put of this edge-detection filter (shown in Figure |2];) marks 
the TRGB. Relatively small photometric errors cause the out- 
put of the edge-detection filter to be noisy. In order to ob- 
tain a smoother output which would yield a more reliable 
distance estimate, we calculate the same edge-detection func- 
tion using photometric errors three times those estimated from 
the artificial star tests. The result is shown in the bottom 
panel of Figure |2] Apparent /' magnitude of TRGB is esti- 
mated to be i'fj^ug = 24.35 ±0.30. Uncertainties are deter- 
mined as half-width-half-maximum of the peak in Figure |2}l. 
With the absolute /'-band magnitude Mi'jrcb = -3.44± 0.10 
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Fig. 2. — Determination of the tip of the red giant branch. I'-band lumi- 
nosity function (a) and smoothed luminosity function (b) are shown in top 
panels before and after completeness correction (dashed and solid lines re- 
spectively). Outputs of the edge-detection filter are shown in the bottom two 
panels - calculated using the photometric error as determined from artificial 
star tests (c). and calculated using a photometric error four times larger than 
calculated from artificial star tests (d); line types are the same as in the top 
panels. 



(lBellazzinill2008 l), we calculate the NGC 7793 distance mod- 
ulus to be m-M = 27.79 ± 0.32, corresponding to the distance 
of 3 .6 1 ± 0. 53 Mpc. This is in good agr eement with the values 
obtained by iKarachentsev et al.' ('2003'. 3.91 ±0.41 Mpc) an d 
the GHOSTS survey (3.73 Mpc Radbur n-Smith et al.ll20lTl) . 

3.2. Color-Magnitude Diagram 

Color-magnitude diagrams (CMDs) of two NGC 7793 are 
shown in Figure [3] Also marked in the Figure are 50% com- 
pleteness limit and photometric errors as determined from ar- 
tificial star tests at a given /' magnitude and color of g'-i' = I. 
The ig' - i',i') CMDs reveal a prominent red giant branch 
(RGB). RGB stars make up more than 60% of detections in 
each field. Region just blueward of the RGB is most likely 
populated by unresolved faint background galaxies. Also po- 
tentially visible in the CMD is a main sequence (MS) and an 
asymptotic giant branch (AGB) population. Selection boxes 
used for deriving properties of individual stellar populations 
are shown in Figure |4] AGB and MS populations make up 

2% and ~ 10% of all objects in the field. 
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Fig. 3. — Color-magnitude diagram of two NGC 7793 fields. Dashed lines 
represent the 50% completeness limit. Photometric uncertainties at a given 
magnitude and color of g' — i' = \ are shown on the left of the panels. 

3.3. Star Counts Profiles 

We derive radial profiles of NGC 7793 in order to examine 
the extent and structure of its outer disk. Star counts pro- 
files are computed by counting up the stars in elliptical annuli 
with inclination and position angle corresponding to that of 
NGC 7793. These raw counts are then normalized by the area 
sampled by each annulus. 

The most challenging aspect of studying outer disk light 
profiles is the need to determine a background level that is 
then subtracted from the raw profile in order to calculate the 
true star counts or surface brightness profile of a galaxy. In the 
case of surface photometry this is manifested through difficul- 
ties in deriving accurate estimates of sky brightness. When 
resolved stellar photometry is used instead to study outskirts 
of spirals, the problem translates into how to reliably evaluate 
the number of unresolved faint background galaxies that are 
mistakenly included into stellar catalogs. Unlike in the case of 
surface photometry, in star counts studies we are not limited 
by the fundamental limits of observations but by the accu- 
racy with which we can determine the galaxy number counts. 
However, the challenge is particularly difficult when informa- 
tion on number counts as a function of color (and not only 
magnitude) is required. We therefore pay careful attention to 
estimating the number counts of faint background galaxies. 

In Figure |5^,e we show star count profiles constructed us- 
ing all stars in the field, irres pective of their colors. >ye use 
the GalaxyCount online tool JeIHs & Bland-Hawthornll2007^ 
to determine galaxy number counts in the magnitude range 
probed by our data; dashed horizontal line in top panels shows 
the estimated background galaxy level of 120(±12) arcmin"^ 
(SE) and 117(±12) arcmin'^ (NW field). The open and full 
symbols in the figure represent the star counts before and af- 
ter the number counts of galaxy contaminants have been sub- 
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Fig. 4. — Color-magnitude diagram of botli NGC 7793 fields with selection 
boxes for RGB, AGB, and MS populations marked. Dashed lines represent 
the 50% completeness limit for SE field (dashed line) and NW field (dotted 
line). 

tracted. Figure [5^ shows that the star counts profile of the 
NGC 7793 NW field falls off to practically zero within the 
extent of our data, at the radius of approximately 11' or 11 .5 
kpc. It is unclear whether the upturn reflected in the two out- 
ermost datapoints for the SE field (Figure|5t) is real, or (pos- 
sibly more likely) a result of a slight underestimation of back- 
ground counts. 

We also calculate the star counts profiles of distinct stel- 
lar populations in the outer disk of NGC 7793 (Figure 
d,f-h). Selection boxes used to separate stars into differ- 
ent populations are marked on the color-magnitude diagram 
in Figure |4] Since the GalaxyCount tool does not pro- 
vide direct information on the color of background galax- 
ies, we require an alternative method for determi ning the 
galaxy number counts. Following lVlaiic et al.l ( 120091) . we em- 
ploy the data from the William Herschel Deep Field (WHDF, 
[Metcalfe et al. 2001) to estimate the contamination from the 
faint background galaxy population. We determine the num- 
ber counts within the asymptotic giant branch (AGB) box di- 
rectly from the WHDF data, while for the stellar populations 
reaching fainter magnitudes than probed with WHDF (i.e. 
main sequence (MS) and re d giant branch (R GB) stars) we 
use the method described in IVlajic et al.l (|2009) to calculate 
the galaxy number counts. We calculate the /'-band number 
counts of all galaxies in WHDF and fit linearly the (log of) 
differential number counts in 0.5 mag bins. In order to deter- 
mine the galaxy counts below the limit of the WHDF survey 
we assume that the counts in the bins 2-3 magnitudes below 
the survey Hmit follow the same Unear trend (in the log space) 



as the counts in the brighter bins used in the fit. We finally 
correct the derived galaxy number counts using completeness 
curves of our data. The resulting background galaxy counts 
are 3.9±2.1, 51 ±7 (44±7) and 52±8 (49±7) arcmin'^, 
for the AGB, RGB and MS selection regions, respectively, for 
the SE (NW) field. (Quoted errors are variance, as estimated 
by GalaxyCount.) While contamination-subtracted profiles 
of RGB stars largely confirm the finding from Figures |5h,e 
(with the distinction that the RGB profile for the NW field 
falls off more steeply and RGB stars are only detected out to 
10') we detect no main sequence stars and all objects within 
our MS selection box can be attributed to the contaminating 
background galaxy population (the galaxy number counts for 
the MS selection box are ^2-3 times higher than the derived 
star counts for this color-magnitude region). We detect AGB 
stars out to 8-9' (8.5-9.5 kpc), after which their number 
counts fall bellow the estimated background galaxy level. As 
we show in Vlaiic et al.. (2009), at the high galactic latitudes 
of the Sculptor G roup, contamination from the Milky Way 
st ars is negU gible ("Robin et al.ll2003l:ISharma et alj|201 lb . 

iRadburn- Smith (2011) find a break in the radial profile 
of young and intermediate age stars in the outer disk of 
NGC 7793 (their HST/ACS fields overlap significantly with 
our SE field), with the scale length of a stellar population be- 
ing shorter for younger stars. This is largely consistent with 
the star counts profiles we derive. Due to the higher level of 
contamination in our ground based data we see no MS stars, 
in agreement with the short scale length for this population 
found bv lRadburn-Smithl(l201 li) : similarly, we find AGB stars 
to be more extended than the MS population, with the RGB 
stars having the largest scale length. 

Comparing the R GB profiles of NGC 7793 and NGC 300 
d Vlaiic et al.ll2009i Figure 9) we find that the counts in the 
outermost bins shown in Figure |5]3 are ^ 2 times lower than 
corresponding counts in the most distant bins in the outer 
disk of NGC 300. This is yet another piece of evidence sup- 
porting our earlier findin g of an ex tended exponential disk 
in NGC 300 (Bland-HawthornetaLi2005; Vlajic et al. 2009). 
While our CMD of NGC 300 reaches 4 mag below the tip 
of the RGB, compared to only 2.5 mag in NGC 7793, back- 
ground galaxy number counts increase rapidly with magni- 
tude and the counts in the faintest magnitude bins domi- 
nate the total galaxy counts. Our CMDs of NGC 300 and 
NGC 7793 reach same apparent depth (^ 26.5-27 mag) and 
hence experience roughly the same contamination by faint 
background galaxies. The difference in star counts in the out- 
ermost bins therefore does not reflect the difference in galaxy 
number counts but in star counts, and is an additional in- 
dependent confirmation of the extended exponential disk in 
NGC 300 out to at least 10 disk scale lengths. 

3.4. Surface Brightness Profile 

The power of resolved stellar photometry over sur- 
face photometry is most easily recognized if star counts 
are transformed into measurements of effective surface 
brightness and compared with existing surface brightness 
data. It has been shown in a number of works recently 
^and-Hawthorn et al. 2005; Irwin et al. 2005; de Jong et al] 
1^71: iRadburn-Smith et alj|201 II) that this approach allows 
one to reach surface brightnesses 3-4 mag arcsec"^ below 
the limit of surface photometry. 

We divide the data in 0.5' wide annuli and calculate surface 
brightness in each annulus as: 
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Fig. 5. — Starcounts profiles of the SE field (left panels) and NW field (right panels). From top to bottom, the profiles shown are of all stars (a & e), RGB (b & 
f), AGB (c & g) and MS stars (d & h). In each panel the dashed line displays the estimated galaxy counts level, determined using methods described in the text. 
Open and full symbols show the star counts before and after the background galaxy contribution has been subtracted. 

where F' = 168887 is a scale factor equivalent to F" that con- 
verts A^y^/t to arcmin^ and Nggi = 120(±12) arcmin"^ for SE and 
1 17(±12) for NW field (ES- We next perform 100 realiza- 
tions of background subtraction; in each of these we remove a 
random set of Ncom stars from each annulus and calculate sur- 
face brightness using the remaining stars as described above. 
The final surface brightness profile is the mean of a hundred 
realizations. We find the standard deviation of the whole set of 
realizations for a given annulus not to be significant (< 0.03 
mag arcsec"^ in g' and < 0.05 mag arcsec"^ in /') confirm- 
ing that the derivation of surface brightness profile is robust 
against the choice of objects we exclude as contamination. 

The derived g' and /' effective surface brightness profiles 
are shown in Figured The in ner disk data points in the figure 
are from the lCarignanI (| 



Here, Npix is a number of pixels in an annulus, fgi and /]/ 
are radial completeness factors and m, are magnitudes of stars 
within a given annulus. Radial completeness of our data is 
lowest in the innermost annulus (45%) due to crowding, and 
increases to an average of 87% in the outermost disk. To con- 
vert surface brightness to units of mag arcsec"^ we multiply 
the effective flux under the logarithm with the inverse of the 
square of the GMOS pixel size (1 pix= 0.146") which is equal 
to F" = 1/0.146^ = 47. 

In order to take into account the contribution from faint 
background galaxies (see ^3.3t we estimate contamination in 
each annulus as 



^gal^ pix 



fHP) study. Smce, to mmimize spuri- 
ous detections, we used the matched stellar catalog rather than 
individual g' and /' catalogs when deriving surface brightness, 
it is not surprising to find that g' and /' profiles are practi- 
cally identical. The uncertainty in the galaxy number counts 
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TABLE 1 

NGC 7793 DISK SCALE LENGTH. 



Field 


Filter 


scale length ['] 


scale length [kpc] 


Carignan 


By 


l.llitO.Ol 


1.17±0.02 


SE 


g' 


1.12±0.05 


1.18±0.05 


SE 


i' 


1.24±0.09 


1.30±0.10 


NW 


g' 


1.12±0.07 


1.18±0.07 


NW 


i' 


1.05±0.06 


1.10±0.07 



Fig. 6. — Surface brightness profile of NGC 7793. Inner disk data (full squares) are surface biightness measurements from Carignaij 419851) . Empty circles 
(squares) are effective surface brightness measurements in g' (/') band, derived from star counts and corrected for background galaxy contamination. Error bars 
reflect uncertainties in galaxy number counts. Due to different relative depths and completeness, data had to be shifted by different amounts to match the lCarignanI 
(1985) measurements; see text for details. Resolved stellar counts reach more than 4 mag arcsec"^ deeper than surface photometry and significantly increase the 
known extent of NGC 7793. 

mination from broad band photometry. 

In calculating the metallicity gradient, we only use stars 
from the top ^ 1 magnitudes of the RGB selection box in Fig- 
ure|4](i.e. with M\ < -2.5). At fainter magnitudes stars at the 
red/metal-rich side of the RGB fall below the completeness 
limit, and the resulting MDF is artificially skewed towards the 
metal-poor end. We adopt four different sets of stellar evolu- 
tionary tracks from' VandenBerg et al.l (I2006h (with ages of 8 
and 12 Gyr and [a/Fe] of 0.0 and 0.3) and interpolate between 
them to derive metallicities on a star by star basis. The model 
grid consist of 16 finely spaced red giant tracks covering the 
range of metallicities from [Fe/H]= -2.31 to 0.00 in the steps 
of approximately 0.1 dex. 

We divide the data into bins of at least 200 stars and for each 
isochrone set determine median metallicities in these bins. 
The metallicity gradient is calculated as a slope of a linear 
fit to the binned metallicities. Effects of using stellar tracks 
with different ages and/or a-enhancements are shown in Fig- 
ure |2] and Table |2] The use of older and more a-enhanced 
isochrones results in lower overall metallicities, but has prac- 
tically no effect on the derived abundance gradient. 

Figure |7] reveals two surprising characteristics - a strong 
positive metallicity gradient in the NW field and overall low 
metallicities in both fields. It becomes clear upon closer in- 
spection of the Figure |7]that the outermost binned datapoint is 
primarily responsible for the overall steep gradient in the NW 
field. The most distant radial bin covers a radial range of ^ 3'; 
a substantial fraction of detections in these outermost few ar- 
cminutes can be attributed to background contamination (Fig- 
ure |5f), making the number count of actual RGB stars very 
low. If the most distant data point is discarded, the abundance 
gradient becomes flat, with the slope and overall metallicity 
consistent with that derived for the SE field. However, the 
value of the slope does not change with the change of the bin 
size, suggesting that the positive abundance gradient in the 
NW field is real. 

On the issue of very low metallicities, we see two possible 
causes for such behavior Errors in the photometric calibra- 
tion would manifest as a displacement of the RGB within the 
color-magnitude diagram, which would in turn result in false 
metallicity estimates. We reject this possibility, since due 
to the initial problems with the photometry (©, photomet- 
ric calibration was checked extensively and independently by 



is reflected in the error bars in Figure |6l Low/high error bar 
limits represent surface brightness profiles calculated using 
background galaxy number counts that are Icr lower/higher 
than number counts used to derive the original profile. In or- 
der to match the inner disk result, effective surface brightness 
data points of the NW field have been shifted downward by 
1.3 and 0.3 mag in g' and /', respectively. For the SE field, 
corresponding values are 0.9 and 0.0 in g' and /'. 

Fitting an exponential to a light profile in Fig ure [6l we 
calculate the scale length of the NGC 7793 disk. ICarignanI 
(|l985) data points yield for the disk scale length a value of 
1 . 11' ± 0.02 (1 . 17 ± 0.02 kpc). The scale lengths derived us- 
ing our data, excluding data points at 10.5 and 11', are pre- 
sented in Table [T] and are largely in agreement with the Carig- 
nan's value. 

In summary, our effective surface brightness derived from 
star counts reaches ~3jnag arcsec"^ deeper than the surface 
brightness profile of lCari gnan ( 1985). We trace the exponen- 
tial light profile out to ^ 9 disk scale lengths, greatly increas- 
ing the known extent of the disk from the Carignan study. 
There is no detectable break in the surface brightness profile 
and the disk remains exponential out to 11' (1 1 .5 kpc). 

3.5. Metallicity Distribution Function and Metallicity 
Gradient in the Outer Disk 

Using the position of stars on the RGB as a proxy for their 
metallicities, we derive the metallicity distribution function 
(MDF) and the metallicity gradient of the stars in the outer 
disk of NGC 7793. Due to higher sensi tivity of the color 
of RG B stars on metallicity than age (e.g. IVandenBerg et al.l 
120061) . the age-metallicity degeneracy does not have a signif- 
icant effect on the shape of the derived MDF and metallicity 
gradient, and allows for a relatively accurate abundance deter- 
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Fig. 7. — M etallicity gradient in NGC 7793. Filled triangle symbols are inner disk [O/H] abundances from the works of lEdmunds & Page]| <1984l) and 
IWebster & Sm ith 1 1983); the solid line shows a linear fit to these abundances. Outer disk metallicity profiles are linear fits to radially binned metallicities 
presented in this work. Different line styles denote abundance profiles derived using different isochi'one sets; ages and [o/Fe] ratios of these isochrones are shown 
in the top right. For clarity, we only show binned metallicities that were used to derive the abundance profile for one isochrone set (filled circles). EiTor bars 
represent mean value of abundance uncertainties in each bin. 



Gemini personnel, and their analysis confirmed our zero point 
values. Alternatively, it is possible that the component we are 
seeing in the outskirts of NGC 7793 is not an outer disk, but a 
galactic halo. However, the effective surface brightness pro- 
file we derive follows the inner disk profile (Figure |6]l and it 
seems very unlikely that the transition from the disk to the 
halo occurs already at 6', wh ere we derive stellar metall icities 
of [Fe/H]w -1 .5. Moreover. iRadburn-Smith et al.1 (|201 1!) de- 
rive metallicities in the outer disk of NGC 7793 and find low 
mean abundances ranging from -1 .23 to -1 .64 in their Fields 
01 and 02 (which largely overlap with our SE field). 

In order to assess the magnitude of the apparent discontinu- 
ity in the region where the inner and the outer disk metallici- 
ties overlap (^ 6 arcmin, Figure|7]l it is necessary to take into 
account (i) the [a/Fe] ratio, which corrects for the fact that in- 
ner and outer disk abundances refer to different elements, and 
(ii) the age-metallicity relation, which reconciles abundances 
in young HII regions with metallicities of the stars on the red 
giant branch which we have assumed to be at least 8 Gyr old. 
Despite numerous results pointing to a flat a ge-metallicity re- 
lationship in the Solar neighborhood (e.g. Edvardsson et alj 
119931) . these works find that the relation steepens for the oldest 
stars. Assuming that the difference in metallicities of young 
HII regions and 8-12 Gyr old RGB stars is - 0.3-0.5 dex 
and [a/Fe]= 0.0-0.3, it is possible to account for the abun- 
dance difference of 0.5-0.8 dex in the transition region in 
Figure|7] At the high end (i.e. 0.8 dex), this makes high inner 
disk metallicities consistent with the very low abundances we 
derive. 

Furthermore, in deriving metallicities of H II reg ions, 
lEdmunds & PageJ dp) and IWebster & Sriiidil Sl98% use 
strong-line abundance indicators. These have been shown 
to result in metallicities which are a factor of two (0.3 
dex) higher than nebular {Tg based) or stellar abundances 
(jBresolin et al. 2009a) and likely enhance the difference be- 
tween the inner (strong-line) and outer disk (stellar) abun- 
dances in Figure |7] 

Similar abundance behavior (although of a smaller magni- 
tude) is observed in the Milky Way. In Figure |8] we repro- 
duce Figure 4 from Carney et al. (2005). Old op en clu sters 
are sh own as diamonds (from Twarog et al.ll 1 997t IChen et alJ 
I2OOI and squares (fr om.Yong et al., 2005^ Circles are three 
field red giants from ICarnev et al.l (120051) . Cepheids from 



TABLE 2 

Metallicity gradient in NGC 7793 . 



Field 


Age (Gyr) 


[a/Fe] 


a 




h 




SE 


8 


0.0 


0.004 


0.011 


-1.49 


0.10 


SE 


8 


0.3 


0.006 


0.011 


-1.69 


0.11 


SE 


12 


0.0 


0.005 


0.012 


-1.64 


0.11 


SE 


12 


0.3 


0.005 


0.011 


-1.82 


0.11 


NW 


8 


0.0 


0.095 


0.023 


-2.17 


0.20 


NW 


8 


0.3 


0.097 


0.024 


-2.37 


0.21 


NW 


12 


0.0 


0.100 


0.023 


-2.35 


0.21 


NW 


12 


0.3 


0.100 


0.024 


-2.53 


0.21 



Note. — Gradients are in the form of [a/Fe]= aR+b, where a 
and the corresponding error have units of dex kpc"' . Abundance 
gradient for the inner disk using the data on [O/H] abundances of 
HII regions from iEdmunds & Paget (il984) and i Webster & SmitH 
fT985) is -O.ll ± 0.02 dex kpc"' 



the wo r ks of lAndrievsky et al.l (|2002bl). lAndrievsky et alj 



^sJcyetai.! (/uuza). lAnarievsKy et aij 

('2002a'), 'Andrievsk vetalJ ( l2002cl) . iLucketalJ (120031) . and 



Andrievsky et al. (2004|) are shown as crosses. Cepheids fol- 
low the well-defined negative gradient out to ^ 15 kpc. On 
the contrary, old stars exhibit steeper gradient in the inner disk 
(out to 10- 12 kpc) and a flattening in the outermost parts. In 
the region 10- 15 kpc there is a clear disconnect of up to 0.4 
dex between the Cepheid abundances and those of old stars, 
similar to what we observe in NGC 7793. 

4. DISCUSSION & CONCLUSIONS 

We can summarize the results presented thus far as follows: 

(i) The outer disk of NGC 7793 is primarily old with red 
giant branch stars as the dominant stellar population and a 
small contribution from asymptotic red giant branch stars. 

(ii) After the contamination from faint background galaxies 
has been taken into account, number counts of RGB stars, as 
well as those derived using a matched catalog of all stars in 
the field, extend out to 10-11', or 10.5-11.5 kpc. 

(iii) The effective surface brightness profile derived from 
star counts traces the disk of NGC 7793 out to 9 disk scale 
lengths and is ^ 3 mag arcsec"^ deeper than the surface pho- 
tometry data of lCarignaiij ( 119851) . Any potential break in the 
light profile may be associated with a specific stellar popula- 
tion, but we see no evidence of a truncation in old stars. 

(iv) The metallicity gradient in the outer disk of NGC 7793 
does not exhibit a negative profile, characteristic of the in- 
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Fig. 8. — Metallicity gradient in tlie disk of tlie Milky W ay. Plott ed are 
data for (i) old open clusters Twarog et al. (from 1997); Ch en et all (from 

f)3) as diamonds and Yong et al. (2005) as squares, (ii) field red giants 
irneyetal. 2005) as circles, and (iii) Cepheids from Andrievsky et al. 
l02a) . ,Andrievskv etaL ( ,2002q) . ,Luck et al., 12003), and Andrievsky et alj 
120041) as crosses. The discrepancy is obvious between old star abundances 
(diamonds, squares and crossed) and Cepheids (crosses) in the region 10— 15 
kpc. 



ner galactic regions. The value of the slope is independent 
of the exact set of stellar evolutionary tracks used (although 
the older, more a-enhanced isochrones result in lower overall 
metallicities). 

(v) The outer disk metallicity gradient is in disagreement 
with the inner disk slope. The inner and outer disk abun- 
dances in the overlap region are potentially in agreement 
after [a/Fe] ratio, age-metallicity relationship and the use 
of specific abundance indicators have been taken into account. 



Results presented in ^3.51 suggest that the abundance gra- 
dient derived from outer disk stellar [Fe/H] metallicities dif- 
fers in slope from the gradient calculated using [O/H] abun- 
dances from inner disk HII regions (Figure [TJ- However, the 
two cannot be directly compared as the latter dataset probes 
recent gas abundances, while the former refers to chemical 
composition in stars that are at least a couple of Gyr old. In 
addition, given the results of Roskar et al. ( 2008a), which sug- 
gest that gas and stellar metallicities are decoupled and follow 
opposite trends, it is difficult to conclude whether our results 
point to an overall abundance gradient that gets shallower or 
steeper with time. (Models of galactic chemical evolution 
are successful in reproduc i ng bo t h trends, e.g. [M olla et al. 
1997'; 'Boissier & Prantzos' '1999'; i Portinari & Chiosii 1 1999c 
Tosi 1988; Chiappini et al. 2001) However , our re sults sup- 
port the scenario presented by Ros kar et al.l (I2008al) in which 
gas abundances become steeper with time (this is consistent 
with relatively steep inner disk slope in NGC 7793) and the 
stellar abundance gradient in old stars is shallower than that 
in young stars. 

Although the stellar metallicity is the primary factor influ- 
encing colors of RGB stars, age-metallicity degeneracy and 
the assumption of single age in calculating the metallicity 
distribution function introduce uncertainties in the derived 
MDF As shown in Figure |7] the derived metallicity gradi- 
ent is practically independent on adopted isochrones, assum- 
ing that age gradient over the extent of the disk is close to 
constant; a non-zero age gradient would result in a different 
metallicity profile. Negative age gradient in the outer disk 
of NGC 7793 would suggest that the real abundance gradi- 
ent has a higher slope than derived under a constant age as- 
sumption. Stellar ages which decrease with radius are indeed 



consistent with the inside-out scenario for galaxy for mation 
dLarsonl 119761: [M atteucci & Francois 1989; Chiappi ni et all 
119971; iNaab & Ostrikerii2006c iMunoz-Mateos et al..i2007h . In 
this picture, a galaxy's inner regions are built up at ear- 
lier times than outer parts, and as a result contain on aver- 
age older stars than outermost regions. However, recent re- 
sults from r esolved stars (Bar ker et al.l 120071; IWilhams et al.l 
120091120 101) and surface photometrv (iBakos et al.l2008l) seem 
to suggest that positive age gradients are frequently ob- 
served in outer disks of spirals. This is supported by re- 
cent simulations of disk evolution (Roskar et al. 2008b al 
iSanchez-Blaizquez et al.li2009l) . which find that radial migra- 
tions of stars within the disk are responsible for the reversed 
age profile at large radii. If the same holds in NGC 7793, 
the true abundance gradient would be negative, flat or mildly 
positive, depending on the magnitude of this effect. 

There is a broad agreement that negative stellar abun- 
dance gradients, easily explained in the context of inside- 
out models for galaxy forrn ation dGoetz & KoeppenI 119921; 
iMatteucci & Franc oisI 49 891). are a common feature of 
disk galaxies (Zaritsky et al.l Il994t [Ferguson et al.l 119981; 
'Gogar ten et al.ll2()10i) . Surface density, yield and star forma- 
tion all decrease with radius, resulting in metallicity distribu- 
tion that is more metal-rich in central parts and decreases pro- 
gressively towards the outer disk. However, abundance pro- 
files in faint outer disks are more difficult to derive and there 
is no general consensus on their shape and origin. Growing 
body of evidence suggests that (most) spirals exhibit a flat- 
tening of their metallicity gradient in the outermost disk. Ob- 
servationally, the strongest case has been made for the Galaxy 
( Andrievskv et al. 2004; Yong et al. 2006; Carraro et al. 2007] 
Pedicefli et al. 2009), M83 (Bresolin et al. 2 009b.). and M3 1 
(Worthev et al. 2005). In the models of Rosk ar et al.1 (l2008b|) . 
IRoskar et all iQQOSa) and ISanchez -Blazqu ez et all (12009^ . 
stellar radial mixing has been shown to be able to produce flat 
abundance profiles by 'smoothing out' the underlying nega- 
tive gradient. On the other hand, a mildl y positive metallici ty 
gradient has been observed in NGC 300 (IVlaiic et al.ll2009h . 

As mentioned earlier, a positive age gradient in the outer 
disk of NGC 7793 would bias our derivation of abundance 
profile and a flat underlying metallicity profile would be ob- 
served as a positive gradient instead. Our positive metallic- 
ity gradient in the NW field could therefore be interpreted as 
a combination of a flat abundance and positive age gradient. 
This particular combination of age and metallicity behavior 
has bee n found to arise as a consequence of stellar migrations 
( IRosk ar et al. 2008a; Sanchez-Blazquez et al. 2009). On the 
other hand, it is possible that the positive metallicity gradient 
in NGC 7793 is real an d does not refle c t the e ffects of age- 
metallicity degeneracy. iMinchev et al.l ( 120101) find that the 
overlap of spiral and bar resonances in the disk triggers sig- 
nificant migration of stars and results in positive abundance 
profile in the outermost regions, similar to what we observe 
in NGC 7793. 

Alternatively, an external mechanism could be responsi- 
ble for the shape of the metallicity gradient in outer disks 
of spirals. In NGC 7793 in particular, the origin of a par- 
ticular abundance profile could be explained by the fact that 
the galaxy harbors a surprisingly small H I disk. While a 
great majority of spirals have more or less extended H I disks, 
sometimes stretching out far beyond the known optical edges, 
neutral hydrogen in NGC 7793 is detected on ly out to 1 1 .5' 
(iCarignan & PuchelfT990t [Walter et al.ll2008h . covering prac- 
tically the same radial extent as our stellar photometry. (In 
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Fig. 9. — (a) Logarithm of surface density in the outer disk of NGC 7793. 
Filled circles are data from Puche et al. 1 1990); solid line is a linear fit to 
these data points, (b) Same as (a), but for rotational velocity, (c) Epicyclic 
frequency derived from the velocity curve in (b). (d) Surface density, (e) 
Toomre Q parameter. 

addition, NGC 7793 exhibits a decreasing velocity curve in 
its outermost parts, which is highly unusual for a galaxy of its 
size.) The reason for a relatively modest H I disk in NGC 7793 
is unclear, particularly given that the galaxy has no obvious 
interactions that could have potentially stripped the gas and 
truncated its H 1 distribution. Evidence for stripped stars in 
the outer disk of NGC 7793 is also lacking. However, it is 
possible to imagine that the upturn in the abundance gradi- 
ent in the NW field is a consequence of a dispersed stream 
of stars that have long fallen below the detectability threshold 
in surface brightness, but still pollute the outer disk metallic- 
ities. Similar to stellar age, positive gradient in [a/Fe] would 
- under the assumption of constant [a/Fe] - be disguised as 
a positive gradient in metallicity. Accounting for the [a/Fe] 
increase from zero to 0.3 over the range covered by our data 
would likely not result in a flat abundance gradient, but would 
certainly lower the slope of the gradient we derive for the NW 
field. 

Finally, the shape of the metallicity gradient in the outer 
disk could be primordial, originating in specifi c galaxy for- 
mation processes taking place at high redshift. iCresci et al] 
( 120101) show that rotating systems at z w 3 show positive abun- 
dance gradients. These are presumed to be generated by cold 
streams depositing pristine material into the centers of galax- 
ies. In today's galaxies the early positive gradient in galac- 
tic center is reversed through processes of star formation and 



chemical evolution, while the signs of the early gradient re- 
main in the outskirts of galaxies. 

From the perspective of potential for star formation or radial 
migration in a given region of a spiral disk, it is interesting to 
examine the radial behavior of the Toomre Q parameter. In a 
thin differentially rotating disk, rotation and pressure work to 
stabilize the disk against axisymmetric perturbations. On the 
other hand, the disk is destabilized by its own surface density. 
The disk is considered unstable against axisymmetric modes 
if Q, given as: 



is less than unity. Here, a is radial velocity dispersion in km 
s"', K is epicyclic frequency in the units of km s"' kpc"', and 
E is gas su rface density in M p, pc~^. We employ H I obser- 
vations of Puche et aTl (119901) to estimate K(r) and S(r). We 
approximate velocity curve and log of H 1 surface density be- 
yond 4' as linear functions of radius (Figure |9^,b) and use 
these fits to calculate k and Q. Epicyclic frequency k is cal- 
culated as: 




and its radial distribution is shown in Figure |9j;. We estimate 
velocity dispersio n in the outer parts of NGC 7793 from the 
THINGS survey (I Walter et all 120081) and adopt a(r) = 6 km 
s"' as a mean velocity dispersion in our observed field. Ra- 
dial distribution of the Q parameter is shown in Figure |9^. 
While using a declining velocity dispersion (rather than con- 
stant as we do here) would result in a slower increase of Q, 
Toomre parameter in the outer disk of NGC 7793 seems to 
be significantly above unity and as such does not support an 
environment in which there is an ongoing large scale star for- 
mation. It is possible however that the azimuthally averaged 
Q profile does not adequately describe star formation in outer 
disk. At very low gas densities in outskirts of spirals, star 
formation might have proceeded in a stochastic manner and 
in small clumps, rather than on global scales addressed by Q 
as calculated here. This is supported by GALEX results of 
clumpy star f ormation in outer disks of nearby spirals (e.g. 
iThilker et al.ll200 7). The high value of Q presents a poten- 
tial problem for radial migration as well. External perturba- 
tions, e.g. from a passing satellite galaxy, have been shown 
to be able to cause r adial mixing in a Milky Way-type galaxy 
(lOuillen et al.ll2009h and could be responsible for stellar mi- 
grations in high-2 environments. 
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